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Abstract 

Transformational optics are shown to markedly enhance the control of the electro- 
magnetic wave trajectories within metamaterials with unconventional functional- 
ities such as a beam splitter, a toroidal carpet, a Luneburg lens and a black hole, 
all of which are specially designed for surface plasmon polaritons propagating on 
a metal plate. 
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1. Introduction 

In 1998, a team led by Ebbesen discovered that resonant excitations of sur- 
face plasmons enhance electric fields at a surface that force light through its tiny 
holes, giving very high transmission coefficients in the sub-wavelength regime 
JT). Five years afterwards, Pendry, Martin-Moreno and Garcia- Vidal subsequently 
proposed a homogenized model of such structured metal surfaces in order to push 
the dispersion relation of surface plasmon polaritons to new boundaries ||2). Plas- 
monics actually holds new promises with recent proposals in invisibility cloaks 
relying upon plasmonic metamaterials that have already led to fascinating results 
01 SI [Hill. These include resonant shells with a suitable out-of-phase polarizabil- 
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ity in order to compensate the scattering from the knowledge of the electromag- 
netic parameters of the object to hide, and external cloaking, whereby a plasmonic 
resonance cancels the external field at the location of a set of electric dipoles. Re- 
cently, Baumeier et al. have demonstrated theoretically and experimentally that 
it is possible to reduce significantly the scattering of an object by a surface plas- 
mon polariton, when it is surrounded by two concentric rings of point scatterers 
(6). In this invited paper, we propose three plasmonic devices designed using 
powerful tools of transformational optics. We validate our theoretical proposals 
with three-dimensional computations using the finite element package COMSOL 
MULTIPHYSICS. This emerging area of photonics is fueled by analogies with 
cosmologic physics in non-Euclidean space-time metrics Q: here, we focus our 
attention on a beam splitter, a toroidal carpet, a Luneburg lens and an optical black 
hole, for surface plasmons propagating at an anisotropic metal-dielectric inter- 
face. These four electromagnetic paradigms emphasize an unprecedented control 
of surface waves using mathematical tools of general relativity (8] |9l in order to 
design some meta- surfaces |fT0l - |fT5l achieving new plasmonic functionalities. 



2. Beamsplitter 

Recent advances in transformational optics have led to the proposal by de 
Rham et al. of a beam splitter which is a heterogeneous anisotropic two-dimensional 
slab splitting an incident beam into two beams propagating along different direc- 
tions. Importantly, the slab is impedanced matched to the surrounding medium, 
hence it does not exhibit any reflection at the entrance or exit boundaries [16]. We 
would like to adapt this device to the case of surface plasmon polaritons propa- 
gating at a transformed metal-dielectric interface. The bottom-neck of the design 
is the following transform lfT6ll : 
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where J rr i is the Jacobian matrix of the transformation, and c\2 is given by : 
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The inverse transformation matrix takes the following form: 
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Figure 1 : Beam splitter, a) Schematic diagram for a SPP beam splitter. A SPP is 
propagating from the left side, b) Top view of beam splitter with regions I and II 

deduced from the geometric transform [T] with y\ = x, y2 = 1 in I and y\ = —x, 
y2 = — 1 in II; c) Three-dimensional plot of the real part of the magnetic field for 
a y-polarised SPP propagating in ^-direction (phase); d) Two-dimensional plot of 

the normalised powerflow of the magnetic field in the (x,y) plane. Stream lines 

appear in red color. 



T_1 =J^J^ det (J^) ( 3 ) 
where J~J denotes the inverse transpose matrix of J xx /, and det(J x y) its determi- 



nent. Besides, the explicit expression of T 1 is: 

(a -1 -c 2 i 
-C2i a(l + c 2 2l ) | (4) 
a 1 

where a = . We show in figure 1 a typical computation for a beam splitter 



in the particular case when the functions y\ and y2 bounding the regions I and II 
are straight lines symmetric with respect to y = (here, y\ = ±x and yi = ± 1 . We 
however note that y\ and y2 could any smooth functions. 



3. Toroidal Carpet 



Next, we propose a design of a toroidal carpet which is deduced from the 
following transform: 

x 1 = x 
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where J zz > is the Jacobian matrix of the transformation as per: 
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and partial derivatives can be expressed as: 
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The resulting transformation matrix is: 
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In the design shown in figure 2, we actually mapped a toroidal bump (de- 
scribed by the surface of altitude zi) on a flat ring (described by the surface zi): 



Zl = ^Jb 2 -( v / x 2 +y 2 -a) 2 +zo with a = 0.3 ; b = 0.1 ; zo = -0.05 
z 2 = 0.2 
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(y / x 2 +y 2 -a)x 
'b 2 - ( a/jc 2 + y 2 - a) 2 . a/jc 2 + y 2 
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(\/x 2 + y 2 -a)y 
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We show in figure 2 the result of our three-dimensional computations, which 
clearly demonstrate that an electromagnetic field can be stored within an invisible 
carpet that does not perturb the ambient SPP field. Applications might be in hard- 
discs for an all-optic computer. 
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Figure 2: Plane wave, a) Electromagnetic field incident from the top (z-axis), 
with the magnitude of the magnetic field being represented; b) Phase 
representation of the magnetic field, y-polarized SPP wave; c) Two-dimensional 
plot in the (x,y) plane; d) Two-dimensional plot in the (x,y) plane; Moreover, an 
electromagnetic field has been launched inside the toroidal carpet filled here with 
a dielectric medium of permittivitty e = 2 and surrounded by a thin metal 
coating. The field circulates in closed trajectories and does not perturb the SPP 

propagation on the metal plate. 
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4. Luneburg lens 



There is currently a renewed interest in gradient index lenses, such as Eaton 
and Luneburg lenses and Maxwell's fisheye. Following the earlier proposal by 
Liu et al. [|T2l of a Luneburg lens which is made of a homogeneous dielectric with 
a specific spatial variation of the profile, we decide here to investigate the case of 
a heterogeneous Luneburg lens, with the usual spatially varying refractive index 
given by lfT8l : 



2-(y / x 2 +y 2 /R) 2 



(8) 



where R is the radius of the lens and the coordinate axis is located at the center of 
the lens. 




Figure 3: Plasmonic cylindrical Luneburg lens with a spatially varying refractive 
index, (a) Schematic diagram for a cylindrical Luneburg lens on a metal plate 
(the color scale ranges from 1 (vacuum), blue color, to 1.414, red color; (b) 
Streamlines show the focussing effect for a SPP incident upon the Luneburg lens 
(view from above); (c) Spatial distribution of the refractive index; (d) Top view 
and side view for a j-polarized SPP propagating from the left and incident upon 
the Luneberg lens at 700 nm (2D plots of the normalized real part of the 
magnetic field in the (x,y) and (x,z) planes). 



We note that the metal is kept untounched, and no magnetism is involved in 
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this design which thus looks as a reasonable design for an experimental vailida- 
tion. 

5. Black hole 
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Figure 4: Plasmonic cylindrical black hole mimicking a Schwarschild metric, (a) 
Schematic diagram for a transformed metal plate with a cylindrical transformed 

medium on the top; (b,c,d) A SPP Gaussian beam making an angle of 60 degrees 
with the x-axis is incident from the right (b: top view, normalized magnitude of 
magnetic field; c: 3D plot for the normalized real part of magnetic field; d: top 

view phase of magnetic field). 



The infinitesimal line element in a Schwarschild metric for a black hole can 
be expressed as (H 

ds 1 = (l - ^)c 2 dtA - (l - dr 2 -r 2 d6 2 -dz 2 , (9) 

where L denotes the event horizon [7], t is the time variable, (r, 0,z) denote the 
space coordinates and c is the velocity of light. Without loss of generality, one can 
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simply skip the time variable within the framework of time-harmonic Maxwell's 
equations. 

Following the proposal of an optical black hole by Chen et al. [19], we red- 
erived the transformation matrix associated with the spatial part of the Schwarschild 
metric and we found that 
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(10) 



where A = 1 — i in the absorbing core (r < L) and A = 1 in the transformation 
layer (r > L). 

In the present case, we launch a SPP at the wavelength X = 700 nm and width 
1400 nm. The radius of the inner absorbing core is L = 0.9 \im and L\ = 1 \im. In 
figure 4, the trajectory of the Gaussian beam is clearly bent around the center part 
of the black hole (see also the streamlines), as it should in a Schwarschild metric. 



6. Conclusion 

In this paper, we reviewed some techniques of transformational optics applied 
to surface plasmon polaritons propagating at a flat interface separating a trans- 
formed Drude metal and a transformed dielectric medium. We have proposed 
four original meta-surfaces: A beam splitter which splits a large Gaussian SPP in 
two narrower Gaussian SPPs with different directions; A toroidal carpet which de- 
tours an SPP without disturbing its wavefront and its amplitude; A Luneburg lens 
which focusses an incident SPP onto a small region behind the lens; An optical 
black hole which can be modeled as vacuum solutions to the Einstein field equa- 
tions. The latter device can be seen as intrinsic parts of the unbounded version of 
the Schwarzschild metric describing an eternal black hole with no charge and no 
rotation. Indeed, all geodesies of a free-falling particle in the spacetime can be 
continued arbitrarily far into the particle's future or past, unless the geodesic hits 
a gravitational singularity like the one at the center of the black hole's interior. In 
a similar way, we have seen that ray trajectories of a time-harmonic solution to 
the Maxwell's equations in the optical black hole converge towards the center of 
the black hole, in some kind of spiral manner, provided the rays intersect the event 
horizon of the black hole. 
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We discovered while finalising the mansucript that a paper had been submitted 
by Zhang's team on numerical and experimental results for a plasmonic Luneberg 
lens, but we emphasize that our results cross-check their simulations. 
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